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ABSTRACT 

A scalar field method to obtain transverse solutions 
of the vector Laplace and Helmholtz equations in 
spherical coordinates for boundary-value problems 
with azimuthal symmetry is described. Neither 
scalar nor vector potentials are used. Solutions 
are obtained by use of separation of variables 
instead of dyadic Green's functions expanded in 
terms of vector spherical harmonics. Applications 
to the calculations of magnetic fields from steady 
and oscillating localized current distributions are 
presented. 

Keywords: Vector Laplace and Helmholtz Equa- 
tions, Spherical Coordinates, Electromagnetic 
Fields. 



1. INTRODUCTION 

In many physical applications the vector Laplace 
equation V 2 F — or the vector Helmholtz equation 
V 2 F + k 2 F = for spherical coordinates and 
boundary conditions has to be solved. It is the 
case for instance with the electric and magnetic 
fields wherever there is no free charge and current 
distributions. In spherical coordinates, however, the 
Laplacian of the vector V 2 F leads to a set of three 
simultaneous equations, each equation involving all 
three components F r , Fg, F v [1]. This complication 
is well known and general techniques for solving 
these equations have been developed, based on a 
dyadic Green's function involving vector spherical 
harmonics which transforms the boundary condi- 
tions and source densities into the vector solution [2] . 

In a number of cases the vector field must also have 
zero divergence: V • F — 0. This is so for example 
for the magnetic field and in general for transverse 
vector fields. With this transverse condition it can 
be shown that the above vector equations can be 
separated into an equation for F r alone. However, 
this is not true for the other two scalar components 
of the vector field. Nevertheless, it is possible 
to obtain a straightforward solution for Fg in a 
problem with azimuthal symmetry; here we are not 
concerned about those symmetry cases where only 



Ftp is nonzero because a scalar technique of sepa- 
ration of variables is already known to obtain the 
solution [3]. Several problems of physical interest 
can then be treated with the standard mathematical 
methods used to solving scalar field equations, so 
avoiding the dyadic method [4] . 

The purpose of this paper is to report the use of 
the technique to compute the solution to instruc- 
tive examples taken from the electromagnetic theory, 
though there are different other places where vector 
fields show up. Specifically, we work out magnetic 
field solutions from localized current distributions. 
The main methods which may be used for solving 
boundary-value problems in magnetism involve the 
vector potential and the scalar potential [5]. The 
new technique we discuss here deals with the mag- 
netic field itself, without reference to any kind of 
potential. The approach can also be used to calcu- 
late electric fields from localized charge distributions 
since the equations to be solved are the same. 



2. TRANSVERSE SOLUTIONS OF THE 
VECTOR LAPLACE EQUATION 

In order to calculate the components of the vector 
V 2 F along the axes of the spherical system of coor- 
dinates, the following identity must be used: 

V 2 F = -V x (V x F) + V(V . p). (1) 

By combining with the transverse condition V--F = 
and requiring azimuthal symmetry, we find the equa- 
tion 



(V 2 F), 
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for the radial component of the vector Laplace equa- 
tion V 2 F = 0. To compute Fg is convenient to use 
the transverse condition: 
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where azimuthal symmetry has been demanded. 

Equation © can be solved by the usual technique of 
separation of variables. If a product form for F r is 
assumed, then it can be written 



F r (r,0) = ^P{6). 



(4) 



The following separate equations for u(r) and P{&) 
are obtained: 
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where n(n + 1) is the real separation constant. The 
solution of Eq. (JSJ is 



u(r) = ar n+i + 



(7) 



Equation JfjJ is the Legendre equation of order n and 
the only solution which is single valued, finite and 
continuous over the whole interval corresponds to the 
Legendre polynomial P„(cos 0), n being a positive 
integer. Thus the general solution for F r is 

Fr(r, 0) = J2( a » r " _1 + ~STa ) P »(cos 9). (8) 

The simplest way of treating Eq. © is to use the 
series expansion 

oo , 

F*(r,0) = 5>»(r) -P„(cos0), (9) 



which leads to Eq. clearly, the solutions given 

in Eqs. ©, © and (HUJ satisfy Eq. (JTSJ). 

In the case of the magnetic field, we recall that the 
boundary conditions at a boundary separating two 
regions are 



n ■ (Bi - B 2 ) = 0, 



n x {Hi - H 2 ) = J s , 



(13) 



where Js is the surface current density and the nor- 
mal unit vector h is drawn from the second region 
into the first one. For a linear medium the consti- 
tutive relation B = [iH holds, [i being the constant 
magnetic permeability. To see how the technique 
works out, we consider the simple, common textbook 
example of the magnetic field produced by the rota- 
tion of a sphere of radius a, uniformly charged with a 
charge Q, with constant angular velocity to, which is 
usually solved with the vector potential method [2,5]. 
The surface current density at r = a is 

J s (a, 6, ip) = sin 9 <p, (14) 
47ra 

where the z axis has been chosen along the rotation 
axis. The problem reduces to using the series solu- 
tions |(HJ , © an d l|10[l of the Laplace equation in the 
regions r > a and r < a to conform to the require- 
ment that the magnetic field must be finite at the 
origin, vanish at infinity and satisfy the boundary 
conditions of Eq. (|T3|) at r = a. Using B = fi H, 
where fj, is the vacuum magnetic permeability, the 
resulting magnetic field can be seen to be 



where v n (r) are functions to be determined. By re- 
placing Eqs. © and © into Eq. ©, it is found that 



v n (r) = — r 
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for n > 1 with a a = 0. The coefficients a n and b n 
are to be determined from the boundary conditions. 

On the other hand, we note that the solution ob- 
tained from Eqs. JSJ), © and (|10|) satisfies the angu- 
lar component of the Laplace equation for transverse 
vector fields and azimuthal symmetry: 

-• 1 d 2 1 d 2 F 
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as expected by consistency. This happens because 
Eqs. @ and © imply 
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This solution describes a uniform magnetic field 
inside the sphere and a dipole field outside with 
dipole moment fa — z LuQa 2 /3. 

Another example is that of the magnetic field from 
a circular current loop of radius a lying in the 
x-y plane with its center at the origin and carrying 
a current /. In this case the surface current density 
at r = a is 

Js(a, 9, ip) = — (5(cos 9) (p. (16) 
a 

The radial component of the magnetic field turns out 



to be now 

n— 



-j2n+2 
„2n+3 



r > a 



'2n 



(17) 



,2n+ 



T ; r < a 



while the angular component becomes 
H e (r,6) = -E ( " in2 " + J )!! ^ +1 (- 
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where -P2n+i( cos 0) = (d/d8) P 2n +i {cos 9) is an as- 
sociated Legendre polynomial. 



3. TRANSVERSE SOLUTIONS OF THE 
VECTOR HELMHOLTZ EQUATION 



F r (r,0) = £ 
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Depending on boundary conditions, the spherical 

(1 2) 

Hankel functions h n : instead of the spherical Bessel 
functions j n , n n may be used. For Fg we again as- 
sume 

Fg (r, 6) = Y. w « W M Pn (C0S 9) ' (23) 
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where for w n we now obtain 
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The constants c„ and d n are computed from the 
requirement that the vector field has to satisfy the 
boundary conditions of the given problem. 

On the other hand, we observe that the solution 
given in Eqs. I|22l) . ll'ol) and (|24[1 satisfies the angular 
component of the Hclmholtz equation for transverse 
vector fields and azimuthal symmetry: 



In the case of the Helmholtz equation for transversal 
vector fields the radial equation becomes 

(V 2 F) r + k 2 F r = ±-^(r 2 F r ) 

r z sin t) oO ou 

Equation © still applies for problems possessing 
azimuthal symmetry; it is used to obtain Fg from F r . 



Next, by putting 

F P (r,0) = ^P(0), 
r 

we can separate Eq. (11911 into two equations: 
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for j(r) and Eq. for P(0), where n(n + 1) is the 
separation constant. Equation (|21fl is the spherical 
Bessel equation of order n in the variable kr. There- 
fore, the general solution for F r is 



(V 2 F) e + k 2 Fg = l^(rF e )--'^ + k 2 F ll 
r or z 
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as expected, where Eq. 1)2 H must be used. 

An illustration is provided by the problem of the 
radiation from a vibrating current loop lying in the 
x-y plane. The surface density current at r = a is 



Js(a,8,tp,t) = — £(cos 
a 
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The series solution of the Helmholtz equation for the 
magnetic field, which is finite at the origin, repre- 
sents outgoing waves at infinity and conforms to the 
boundary conditions of Eq. (|13fl at r = a, is then 
given by 



H r (r,9,t) = i- 
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where the upper lines hold for r > a and the lower 
lines for r < a. The radiative part of the magnetic 
field takes on the limiting form 
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The resulting magnetic field can be seen to be the 
same as that obtained by another more extensive 
methods [2]. 



4. CONCLUSION 

The expressions for the components of V 2 F in 
spherical coordinates are really messy but they 
appear for instance in Laplace and Helmholtz equa- 
tions for vector fields, such as the electromagnetic 
fields. We have presented a fairly simple technique 
of separation of variables to compute transverse 
solutions with V • F = of these vector equations 
applied to boundary-value problems with azimuthal 
symmetry, so avoiding the general, more extensive 
method in which dyadic Green's functions and 
series expansions in vector spherical harmonics 



are introduced. It involves scalar solutions for 
the spherical components of F which resemble the 
ones obtained for the scalar Laplace and Helmholtz 
equations. No scalar or vector potentials are used. 
Illustrations were provided by problems concerned 
with the magnetic field from steady and oscillating 
localized currents. We finally remark that the 
method can also be extended to solving boundary- 
value problems in cylindrical coordinates possessing 
azimuthal symmetry, although in this case we deal 
with a separated equation for each scalar component 
of the vector field [1], so that the problem of solving 
the vector Laplace or Helmholtz equation reduces 
to that of solving the corresponding scalar equation. 
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